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Abstract 

Spin relaxation time of conduction electrons through the Elliot- Yafet, 
D'yakonov-Perel and Bir-Aronov-Pikus mechanisms is calculated theoretically 
for bulk GaAs, GaSb, InAs and InSb of both n- and p-type. Relative im- 
portance of each spin relaxation mechanism is compared and the diagrams 
showing the dominant mechanism are constructed as a function of temper- 
ature and impurity concentrations. Our approach is based upon theoretical 
calculation of the momentum relaxation rate and allows understanding of the 
interplay between various factors affecting the spin relaxation over a broad 

range of temperature and impurity concentration. 
PACS numbers: 72.25.Rb, 76.20.+q, 76.60.Es 
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I. INTRODUCTION 



Recently, intensive experimental and theoretical efforts have been concentrated on the 
physics of electron spin due to the enormous potential of spin based devices. In these 
so called "spintronic" devicesJ3~i information is encoded in the spin state of individual 
electrons, transferred with the electrons, and finally put under measurement. Electron spin 
states relax, i.e., depolarize, by scattering with imperfections or elementary excitations such 
as other carriers and phonons. Therefore, to realize any useful spintronic devices, it is 
essential to understand and have control over spin relaxation such that the information is 
not lost before a required operation is completed. 

The investigation of spin relaxation has a long history dating back to the fifties and 
most studies have concentrated on III-V semiconductors since direct measurement of spin 
relaxation time is possible through optical orientation in these materials. Three main spin 
relaxation mechanisms, the Elliot- Yafetll (EY), D'yakonov-Pereli (DP) and Bir-Aronov- 
PikusS (BAP) mechanism have been suggested and confirmed experimentally. Earlier works 
for spin relaxation are mainly on bulk systems such as p-GaAsJi~0 p-GaSb,0 GaAlAs^l 
and n-InSb.lii More recently, spin relaxation has also been investigated in quantum well 
structures (GaAs,0 GaAsSb,0 InGaAs/InGaAsP0 and GaAs/AlGaAsEl) as well as in bulk 
systems (n- GaAsSS and InAsi). On the theoretical side, there are recent approaches 
which refine or extend the original calculations of Refs. 3 and 4 to explain newly obtained 
experimental results. Flatte and coworkersSil employed a nonperturbative 14-band cal- 
culation for the DP mechanism both for bulk and quantum well structures and achieved a 
better agreement with the experimental results. The BAP process was reconsidered through 
a direct Monte Carlo simulation and extended to quantum wells by Maialle and coworkers.il 

In the most studies, the strategy has been to find the relevant spin relaxation mechanism 
by comparing experimental results for spin relaxation time, t s , with theoretically predicted 
dependence on temperature or doping concentrations. Based upon these results, a "phase 
diagram" -like picture showing the dominant spin relaxation mechanism can be constructed 
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to provide a comprehensive global understanding for competition of spin relaxation mecha- 
nisms. However, since available experimental results for r s are usually limited to a narrow 
range of external physical parameters except some intensively investigated materials, such 
pictures are currently available only for p-GaAs and p-GaSb.0 

In this paper, we calculate the electron spin relaxation time for the EY (t^ y ), the 
DP (t s dp ) and the BAP ij^ AP ) processes for several bulk III-V semiconductors: GaAs, 
GaSb, InAs and InSb of both n- and p-type. Our result for t s is based upon theoretical 
calculation of the momentum relaxation time, t p . A diagram is constructed illustrating the 
dominant spin relaxation processes as a function of temperature and impurity concentration 
for each material. The resulting "phase diagrams" for p-GaAs and p-GaSb are in qualitative 
agreement with that of an earlier study.0 The diagrams for the other materials considered 
in this work were not available in the literature and represent an attempt to provide a better 
understanding of interplay between various factors for t s . We also discuss some incomplete 
aspects of the current theories for spin relaxation. 

The rest of this paper is organized as follows. In Sec. II the basic formulation of the 
three spin relaxation mechanisms is briefly described. The details of our calculation for the 
momentum relaxation time (r p ) and r s are presented in Sec. III. In Sec. IV the results for 
r s are compared with available experimental results and the "phase diagrams" for dominant 
spin relaxation is constructed. The conclusion follows in Sec. IV. 



II. RELEVANT SPIN RELAXATION MECHANISMS 

A. Elliot- Yafet Mechanism 

The EY mechanism originates from the fact that in the presence of spin-orbit coupling, 
the exact Bloch state is not a spin eigenstate but a superposition of them. This induces 
a finite probability for spin flip when the spatial part of electron wavefunction experiences 
a transition through scattering even if the involved interaction is spin independent .Si The 
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spin relaxation time is given by@ 



1 . (k B T\ 2 „ (I -v/2\ 2 1 



where £7 S is the band gap and 77 = A/(E g + A) with the spin-orbit splitting of the valence 
band A. A is a dimensionless constant and varies from 2 to 6 depending on dominant 
scattering mechanism for momentum relaxation. 

B. D'yakonov-Perel Mechanism 

In III-V semiconductors, the degeneracy in the conduction band is lifted for k ^ due 
to the absence of inversion symmetry. The resulting energy difference for electrons with the 
same k but different spin states plays the role of an effective magnetic field and results in 
spin precession with angular velocity u;(k) during the time between collisions. Since the 
magnitude and the direction of k changes in an uncontrolled way due to electron scattering 
with impurities and excitations, this process contributes to spin relaxation. This is called 
the DP mechanism^ and r s DP is given byi'0 

r?r ~ W h 2 E g p ' [ ) 

where Q is a dimensionless factor and ranges 0.8—2.7 depending on the dominant momentum 
relaxation process, a is the parameter characterizing the fc 3 -term for conduction band 
electrons and is approximately given by 

Ar] m c 



a 



(3) 
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Here m c and m are the effective mass of the conduction electron and the electron rest mass, 
respectively. 

C. Bir-Aronov-Pikus Mechanism 



Electron spin flip transition is also possible by electron-hole scattering via exchange and 
annihilation interactions. This is called the BAP mechanism and is especially strong in 
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p-type semiconductors due to high hole concentrations. r^ AP is given by several different 
expressions depending on the given external parameters. In the case of a nondegenerate 
semiconductor^!! (N A < N c ), 



2a| /2e\ 1/2 



rf AP t v b \m c 
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where n a j (n at b) is the concentration of free (bound) holes and N c is the critical hole con- 
centration between degeneracy and nondegeneracy. e is the conduction electron energy and 
To is given by the relation 
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with A exc the exchange splitting of the exciton ground state. a B , v B and E B are defined as 

h 2 e _ fm \ 
e z m R \m R J 

h 

v B = ; 

m R a B 

h 2 (m R \ K 

where m R is the reduced mass of electron and hole, a® the Bohr radius (~ 0.53 A) and 7Z the 
Rydberg constant (~ 13.6 eV). ^(r) represents wavefunction describing the relative motion 
of electron with respect to hole and |-^(0)| 2 is the Sommerfeld factor given by 



\m\ 2 = -(i-e-^T\ «=J4- 

k \ E B 



For degenerate case (Na > N c ), the result isL 

1 2o| / e \ x(2e/m c ) 1 / 2 ; if e f < e(m v /m c ), 



-BAP 



n a \mr{ (5) 



TqVb V e // \ x(2e f /m v ) 1 / 2 ; if ef > e (m v /m c ), 

where m v is the hole effective mass and e/ the hole Fermi energy, 



III. CALCULATION 

We first calculate the momentum relaxation time r p . We include contributions from 
the polar optical phonon scattering (t^°), ionized impurity scattering (rjf), piezoelectric 
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scattering (r^ e ), and acoustic phonon deformation potential scattering (r*). Our calculation 
of r p is performed with three simplifying assumptions: 

(a) the classical Boltzmann statistics is assumed for conduction electrons, 

(b) the electrons are scattered in a parabolic band, 

(c) the Mathiessen's rule is applied so that l/r p = I/7* + l/r™ + I/t^ + 1/r*. 

Under these assumptions, r p can be obtained in a straightforward way for the given material 
parameters of a III-V semiconductor. 

According to the Ehrenreich's variational calculation^ is obtained as 

4 h ( eoe^ \ /moV/ 2 e^/ T -l wl) 



j-PO 



-G<-%-*, (6) 



p 3^VKt^T Veo-eoo/ \mj 6 t /T 
where eo and are the low- and high-frequency dielectric constants. Q\ is the longitudinal 
optical phonon frequency converted in the unit of temperature and G^e~^ is calculated as 
in Ref. 26 as a function of temperature and the free carrier density n. 
Tp is described by the Brooks-Herring equation^ 

u _ 1 el/al K{k B T)W fmA W xe~* 



v 3rr 3 / 2 2N m + n W>/ 2 \mj Jo g(n,T,x) 

where N m is the concentration of minority impurities, i.e., acceptors for n-type and donors 
for p-type, and a; is a dimensionless quantity representing (e/ksT). g(n,T,x) is given by 

g(n, T, x) = ln(l + b) - 6/(1 + b) 

with 
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after spherical average of the piezoelectric and elastic constants over the zinc-blende structure 
is performed.il Here hu is the one independent piezoelectric constant and q and q are the 
average longitudinal and transverse elastic constants given by 



Q = (3cu + 2ci2 + 4c 44 )/5; 
Q = (cu - C12 + 3c 44 )/5. 

Finally, Bardeen and Shockley0 showed that r^ p is given by 

dp _ 8^ Ml 5 / 2 /m y/2 a 3 Q 

p 3 E 2 (k B T) 3 / 2 \m c J K ' 1 J 

where i?i is the deformation potential. 

The free carrier concentration n (i.e., electrons for n-type and holes for p-type) is calcu- 
lated from the equation 

n(n + N m ) _N{T) 



N M -N m -n 2 rK k B T' 

Here, Nm is the majority impurity concentration. N{T) is given by {2mk B T/(irh 2 )] 3 / 2 /4:, 
where m represents m c for n-type and m v for p-type, respectively. E^ is the ionization energy 
for majority impurity and is given by {1Z/ t^)(m/mo). 

Table I shows the values of material parameters used in the calculation of t p and t s . 
E g (T) is obtained by linearly interpolating or extrapolating E 9i i and E 9t h and N m is fixed to 
5 x 10 13 cm -3 in most cases. Figure 1 plots the results of mobility calculation, /i = (e/m c )r p , 
for n-GaAs and n-InAs. Good agreement is obtained with the published result of Rode and 
KnightS for n-GaAs while our result for n-InAs shows a larger discrepancy up to ~50% with 
those of Rodeil. This seems to result from the fact that the nonparabolicity of conduction 
band, which we neglected, is stronger in InAs. 

Figure 2 illustrates the dominant momentum relaxation mechanism for n-GaAs as a 
function of temperature and impurity concentration. It is found that the contribution from 
the polar optical phonon scattering is dominant for the high-T and lightly-doped regime, 
while the ionized impurity scattering dominates otherwise. The same qualitative features 
are found for all other materials investigated, both for n- and p-type. 

As was noted previously, both rf Y and t® p include dimensionless factors, i.e., A in 
Eq. ([!]) and Q in Eq. (Q), which vary depending on the dominant momentum relaxation 
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process. At the current stage, it is not clear how the crossover behavior is given quantitatively 
when there is a switch between two momentum relaxation processes. Therefore, we fix the 
dimensionless constants to their median values, i.e., A = 4 and Q = 1.75. This introduces 
~50 % uncertainty in our result for rf Y and rf p . One might correct this error by directly 
looking into the dominant momentum relaxation process. 

To calculate rf AP , we first need to identify the adequate regime for a given parameter 
set. N c is determined by the Mott criterion^ N c (0.26/a^) 3 where a# = aoeo / (m v / m ) . 
The thermal averaged value of 1/t pap is obtained as 

assuming a classical Boltzmann distribution for conduction electrons. On the other hand, 
the expressions for 1/t^ y and 1/t s £>p in Eqs. ([3]) and (|2]) are after thermal averaging with 
respect to e. A difficulty with the calculation of r PAP lies in the fact that there is no reliable 
data for A exc , on which t pap has the dependence of ~ 1/A^ C , for p-InAs and p-InSb. 
Therefore, we examine the tendency of t pap as a function of A exc as well. 



IV. RESULTS AND DISCUSSION 

We first compare the relative importance of each spin relaxation mechanism. Figure 3 
shows the dominant spin relaxation processes for n-type GaAs, GaSb, InAs and InSb. For 
n-type semiconductors, the contribution of the BAP mechanism is negligible since the equi- 
librium hole concentration is extremely small. Therefore, we watch the competition between 
the EY and the DP processes. As shown in Fig. 3, it turns out that for all materials inves- 
tigated there exists a transition from the DP-dominant regime to the EY-dominant regime 
at T <~5 K as the temperature is lowered. These results are consistent with the previ- 
ously published results that the DP process is the relevant spin relaxation mechanism for 
n-GaAs@il and n- InAs0 at high temperature of 300 K and that the EY process is rele- 
vant for n-InSb at low temperature of T = 1.3 K.0 When the acceptor, i.e., the minority 
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impurity, concentration decreases, we find that the DP-dominant regime enlarges. This can 
be understood from following consideration. The acceptors in n-type materials are always 
ionized and the decrease in the acceptor concentration corresponds to the decrease in the 
number of scattering centers for ionized impurity scattering procedure, which is the main 
momentum relaxation mechanism at low temperature. Therefore, a larger r p results as the 
acceptor concentration decreases and this induces a larger rf Y and a smaller t® p since 
rf Y ~ r p and rf p ~ l/r p . 

The diagrams for p-type materials are illustrated in Fig. 4 with 10 14 cm -3 < Na < 10 20 
cm -3 and Np = 5 x 10 13 cm -3 . For p-type materials, no systematic changes are found when 
the minority carrier concentration is varied. For p-GaAs and p-GaSb, we find that the BAP 
(DP) is dominant in the low-T (high-T) and high (low) doping regime. This is in qualitative 
agreement with the results of Aronov et al, 111 in which diagrams of the same idea were 
constructed based on experimental results. For p-InAs, a similar feature to those of p-GaAs 
and p-GaSb is found for A exc = 10 eV, and as A exc decreases, the BAP dominant regime 
becomes smaller. For p-InSb, we obtain similar results to those for p-InAs as a function of 
A exc . Figure 4(d) shows the case of A exc = 0.2 eV where a BAP-dominant regime exits at 
T < 100 K and intermediate doping concentrations. We find abrupt discontinuities in t? AP 
at Na = N c , which results in unphysical sharp cusps at Na ~ 10 18 cm -3 in Fig. 4. This is 
an artifact resulting from the fact that no quantitative expression for l/r PAP is available 
for the crossover between nondegenerate [Eq. (|j)] and degenerate [Eq. ([5])] hole regimes. 
Experimentally, it was found that there exists an intermediate regime at Na ~ N c where r s 
remains nearly flat with respect to the change in Na and that the range of such intermediate 
regime varies depending on the material.0 

Figures 5 and 6 provide the total spin relaxation time, i.e., r s = (l/rf Y + l/r® p )~ l for 
n-type and r s = (l/r PY + l/r s DP + l/r s BAP ) _1 for p-type, respectively. r s ranges from 1 ps 
to 100 ns for n-type materials and from 0.1 ps to 10 ns for p-type materials, respectively, 
over the parameter space shown in Figs. 5 and 6. For n-type materials, t s increases as T 
decreases with the longest r s found at Np ~ 10 17 — 10 18 cm -3 instead of in purer materials. 



This is because the regime shown in Fig. 5 is dominated solely by the DP-process and 1/t^ p , 
which is proportional to r p , increases as the impurity concentration decreases. The same 
qualitative feature has also been found in a recent experiment .0 In our result for n-GaAs, 
r s ranges from 5 ns to 60 ns for T = 25 K, which gives a reasonable agreement with the 
experimental result of Ref. 19 (r s ~ 70 ns at T = 20 K). As for n-InAs with Np = 10 16 
cm~ 3 and T = 300 K, our result gives t s = 12 ps which compares very well with a recent 
experimental result of t s = 19 ± 4 ps.0 

At lower temperature, we find a discrepancy with recent experimental result for n-GaAs. 
In experimental r s ~ 100 ns at 5 K for No = 10 16 cm" 3 was reported, while our result 
predicts a larger value of t s ~ 6 x 10 3 ns. Reference. 19 suggested that the main spin 
relaxation at this low temperature regime is due to the EY mechanism. According to our 
result, however, since r p ~ 1 ps and rf Y and r s DP are given by 7 x 10 4 ns and 6 x 10 3 ns, 
respectively, neither the EY nor the DP mechanism provides a satisfactory explanation for 
the experimental result. Very recently, in ref. 20, spin relaxation time of 290 ± 30 ns at 4.2 
K was reported for bound electrons to donors in n-GaAs and the relevant spin relaxation 
mechanism was proposed to be the hyperfine interaction with nuclei,^ which was not taken 
into account in our current work. A further research incorporating this effect is needed to 
resolve the discrepancy between our result and the experimental result of Ref. 19. 

In p-type materials, smaller r s , i.e., stronger spin relaxation rate, than that in n-type 
materials is found due to the effect of the BAP process. The strong discontinuities at 
Na = N c are also noticeable in Fig. 6 due to the incompleteness of the BAP expressions 
given by Eqs. and (H), as mentioned earlier. 

V. CONCLUSION 

In this paper, we calculated theoretically t s for several bulk III-V semiconductors and 
compared the contributions from the three main spin relaxation mechanisms as a function 
of temperature and donor /acceptor concentrations. In n-type materials, the DP mechanism 
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is found to be dominant down to very low temperature, below which the EY mechanism 
dominates. While our calculated spin relaxation times are in reasonable agreement with the 
experimental results for high temperature regime of T >~20 K, there exists a discrepancy 
at T ~ 5 K for n-GaAs. Further theoretical efforts incorporating other spin relaxation 
mechanisms neglected in this paper are needed for its resolution. As for p-type materials, 
the BAP (DP) mechanism is dominant at low (high) temperature and high (low) acceptor 
concentrations. We find that the crossover between various regimes for spin relaxation re- 
quires a further theoretical investigation for a more thorough understanding and realistic 
comparison with experimental data. This is especially the case for the crossover between 
nondegenerate and degenerate hole regimes for the BAP process. 

Acknowledgment 

We are thankful to M. I. D'yakonov for critical comment and to J. M. Kikkawa for 
useful discussion. This work is supported by the Office of Naval Research and the Defense 
Advanced Research Projects Agency. 



11 



REFERENCES 

X D. DiVincenzo, Science 270, 255 (1995). 
2 G. Prinz, Phys. Today 48, 58 (1995). 
3 L. Sham, Science 277, 1258 (1997). 
4 R. J. Elliot, Phys. Rev. 96, 266 (1954). 

5 Y. Yafet, in Solid State Physics, edited by F. Seitz and D. Turnbull (Academic, New York, 
1963), Vol. 14. 

6 M. I. D'yakonov and V. I. Perel, Sov. Phys. JETP 33, 1053 (1971); Sov. Phys. Solid State 
13, 3023 (1972). 

7 G. L. Bir, A. G. Aronov, and G. E. Pikus, Sov. Phys. JETP 42, 705 (1976). 

8 A. H. Clark, R. D. Burnham, D. J. Chadi, and R. M. White, Solid State Commun. 20, 
385 (1976). 

9 G. Fishman and G. Lampel, Phys. Rev. B 16, 820 (1977). 

°V. I. Maruschak, M. N. Stepanova, and A. N. Titkov, Sov. Phys. Solid State 25, 2035 
(1983). 

1 K. Zerrouati, F. Fabre, G. Bacquet, J. Bandet, J. Frandon, G. Lampel, and D. Paget, 
Phys. Rev. B 37, 1334 (1988). 

2 A. G. Aronov, G. E. Pikus, and A. N. Titkov, Sov. Phys. JETP 57, 680 (1983). 

3 A. H. Clark, R. D. Burnham, D. J. Chadi, and R. M. White, Phys. Rev. B 12, 5758 
(1975). 

4 J. N. Chazalviel, Phys. Rev. B 11, 1555 (1975). 

5 Y. Ohno, R. Terauchi, T. Adachi, F. Matsukura, and H. Ohno, Phys. Rev. Lett. 83, 4196 
(1999). 

12 



K. C. Hall, S. W. Leonard, H. M. van Driel, A. R. Kost, E. Selvig, and D. H. Chow, Appl. 
Phys. Lett. 75, 4156 (1999). 

J. T. Hyland, G. T. Kennedy, A. Miller, and C. C. Button, Semicond. Sci. Technol. 14, 
215 (1999). 

A. Malinowski, R. S. Britton, T. Grevatt, R. T. Harley, D. A. Ritchie, and M. Y. Simmons, 
Phys. Rev. B 62, 13034 (2000). 

J. M. Kikkawa and D. D. Awschalom, Phys. Rev. Lett. 80, 4313 (1998). 

R. I. Dzhioev, B. P. Zakharchenya, V. L. Korenev, D. Gammon, and D. S. Katzer, JETP 
Lett. 74, 182 (2001). 

T. F. Boggess, J. T. Olesberg, C. Yu, M. E. Flatte, and W. H. Lau, Appl. Phys. Lett. 77, 
1333 (2000). 

W. H. Lau, J. T. Olesberg, and M. E. Flatte, Phys. Rev. B 64, 161301 (2001). 

M. Z. Maialle, Phys. Rev. B 54, 1967 (1996); M. Z. Maialle and M. H. Degani, Appl. 
Phys. Lett. 70, 1864 (1997); M. Z. Maialle and M. H. Degani, Phys. Rev. B 55, 13771 
(1997). 

G. E. Pikus and A. N. Titkov, in Optical Orientation, edited by F. Meier and B. P. 
Zakharchenya (North-Holland, Amsterdam, 1984). 

H. Ehrenreich, Phys. Rev. B 120, 1951 (1960). 

H. Ehrenreich, J. Phys. Chem. Solids 8, 130 (1959). 

H. Brooks, Adv. Electron. Electron Phys. 7, 158 (1955). 

H. J. G. Meijer and D. Polder, Physica 19, 255 (1953). 

J. D. Zook, Phys. Rev. 136, A869 (1964). 

J. Bardeen and W. Shockley, Phys. Rev. 80, 72 (1950). 

13 



P. P. Edwards and M. J. Sienko, Phys. Rev. B 17, 2575 (1978). 
D. L. Rode and S. Knight, Phys. Rev. B 3, 2534 (1971). 

D. L. Rode, Phys. Rev. B 3, 3287 (1971). 

Landolt-Bdrnstein Numerical Data and Functional Relationships in Science and Tech- 
nology, edited by O. Madelung, M. Schultz, and H. Weiss, New series, Vol. 17 (1982); 
reprinted in O. Madelung, editor, Semiconductors- Basic Data, 2nd. ed. (Springer, New 
York, 1996). 

E. Haga and H. Kimura, J. Phys. Soc. Japan 19, 658 (1964). 
D. D. Sell, Surf. Sci. 35, 863 (1973). 

A. Abragam, The Pricnciples of Nuclear Magnetism (Clarendon, Oxford, 1961). 



14 



TABLES 

TABLE I. Material parameters. N c is from the relation N c m (0.26/a^) 3 and all other numbers 
are from Ref. 34 unless specified otherwise. 





GaAs 


GaSb 


InAs 


InSb 


m c /m 


0.065 


0.0412 


0.023 


0.0136 


m v /m 


0.5 


0.28 


0.43 


0.45 


A (eV) 


0.341 


0.75 


0.38 


0.85 


E g ,i (eV) 


1.52 (0 K) 


0.822 (0 K) 


0.418 (4.2 K) 


0.235 (1.8 K) 


E g ,h (eV) 


1.42 (300 K) 


0.75 (300 K) 


0.354 (295 K) 


0.23 (77 K) 


eo 


12.515 


15.69 


15.15 


16.8 


Coo 


10.673 


14.44 


12.25 


15.68 


Oi (K) 


410 


335 


343 


280 


en (dyn/cm 2 ) 


1.221xl0 12 


8.834xl0 n 


8.329xlO n 


6.669xlO n 


C12 (dyn/cm 2 ) 


5.66xl0 n 


4.023xl0 n 


4.526xlO n 


3.645xlO n 


C44 (dyn/cm 2 ) 


5.99xl0 n 


4.322xl0 n 


3.959xlO n 


3.02xlO n 


/114 (V/cm) 


1.45xl0 7 


9.5xl0 6 


3.5xl0 6 


4.7xl0 6 


£1 (eV) 


6.3 a 


6.7 6 


4.9 6 


7.2 6 


A exc (fjieV) 


50 c 


24^ 


unknown 


unknown 


N c (cm" 3 ) 


7.53xl0 18 


6.71xl0 17 


2.7xl0 18 


2.27xl0 18 



Ref. 35. b Ref. 33. c Ref. 36. d Ref. 12. 
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FIG. 1. Mobility vs. temperature for (a) n-GaAs for No = 10 14 cm 3 and Na = 5 x 10 13 cm 3 
and (b) n-InAs for N D = 2 x 10 16 cm" 3 and N A = 5 x 10 13 cm . The lines are our calculation 
and the points are from (a) Rode and KnightH and (b) Rode.il 
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FIG. 2. Dominant momentum relaxation process for n-GaAs as a function of temperature and 
donor concentration with Na = 5 x 10 13 cm~ 3 . is in cm~ 3 . 
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FIG. 3. Dominant spin relaxation mechanism for n-type materials. The higher temperature 
regime is governed by the DP mechanism as shown while the lower temperature regime by the EY 
mechanism. Njj is in cm -3 and Na is fixed to 5 x 10 13 cm' 3 . Material parameters are as specified 
in Table I. 
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FIG. 4. Dominant spin relaxation mechanism for p-type materials. Na is in cm" 3 and Njj is 
fixed to 5 x 10 13 cm -3 . A exc = 1, 3 and 10 [j,eV from bottom to top for p-InAs and is fixed at 0.2 
^eV for p-InSb. Other material parameters including A exc for GaAs and GaSb are as specified in 
Table I. 
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FIG. 5. Total spin relaxation time for n-type materials. Np is in cm and t s is in second. Na 
is fixed to 5 x 10 13 cm -3 . 
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FIG. 6. Total spin relaxation time for p-type materials. Na is in cm~ 3 and t s is in second. Np 
is fixed to 5 x 10 13 cm" 3 and A. exc to 1 /ieV for p-InAs and p-InSb. 



21 



